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We investigate SU(2) lattice gauge theory in four dimensions in the max- 
imally abelian projection. Studying the effects on different lattice sizes we 
show that the deconfinement transition of the fields and the percolation tran- 
P^ \ sition of the monopole currents in the three space dimensions are precisely 

related. To arrive properly at this result the uses of a mathematically sound 
characterization of the occurring networks of monopole currents and of an 
appropriate method of gauge fixing turn out to be crucial. In addition we 
investigate detailed features of the monopole structure in time direction. 



1. Introduction 



For the conjecture |jl|, ||] that condensing magnetic monopoles describe quark con- 
finement nonperturbative support first came from compact U(l) lattice gauge the- 
ory P, 3 . To treat the nonabehan case 't Hooft [Q proposed to extract the relevant 
abelian degrees in SU(N) gauge theory by a suitable gauge fixing procedure called 
abelian projection. A nonperturbative realization of this concept has become possi- 
ble on the lattice [0. Among the possible choices of such projections the maximally 
abelian projection |^ appears most appropriate and has been considered in numer- 
ous papers |§]. 

To confirm this picture in Monte Carlo simulations it is attractive not only to 
measure observables but to analyze the configurations and to look for the structures 
of monopole currents which are characteristic for the phases and thus are able to 
signal the deconfinement transition. First analyses of this type have already been 
presented some time ago [^, |iy]. A consideration of the connection of two points by 
currents [^ is also in this spirit. Analyses of the sizes of the structures have been 



given in [|r^, |TB[. More recently a quantity called "wrapping number" by the author 
ll4| has been used in the deconfinement phase. 



Monopole currents which (on lattices with periodic boundary conditions) wrap 
around the torus have been typically observed in the confining phase. However, as 
soon as one goes into more detail and tries to make the characterization precise, 
the problem becomes obvious that one has to deal with complicated networks of 
monopole currents rather than with simple loops, for which the topological charac- 
terization would be straightforward. Decomposing the networks into loops, apart 
from being highly ambiguous, is not allowed because it changes the topology. 

In the context of the deconfinement transition this problem so far has not been 
addressed. In [|l^ it has been noted that the considerations there could not be 
extended to the confining phase because of the occurrence of entangled structures 
rather than of simple loops. Otherwise the indicated problem so far has not been 
mentioned in the respective literature. 

In U(l) gauge theory in four dimensions the related problem has recently been 
solved using mappings which preserve homotopy |T^, 0. This has lead to an un- 



ambiguous characterization of the phases by the presence or absence of a current 
network which is topologically nontrival in all directions. More generally this means 
presence or absence of an infinite network where the definition of "infinite" on the 



finite lattice is dictated by the boundary conditions IT^- For periodic boundary 



conditions it has been checked in detail that neither size nor extension but precisely 
the topological properties are the crucial ones. 

These results suggest to try an analogous characterization of the deconfinement 
transition. In that case the percolation phenomenon of the monopole currents should 
occur only in the three space directions. Due to the fact that the characterization 
is by the indicated topological properties there is no problem giving a consistent 
prescription also in the three-dimensional subspace. 

Within the latter respect a conceptual problem exists for the other quantities 



considered so far, i.e. for monopole density M, O, sizes of monopole loops ||T2l 13 



and quantities related to cluster size |T0[ . We have checked, using larger lattices and 
higher statistics, that such quantities also empirically do not lead to a convincing 
char act erizat ion . 

The indicated characterization by the topology of networks has the advantage that 
it is insensitive to short range fluctuations. In fact, as has been shown in the case 



of U(l) gauge theory [T^, |T6|, ^7j, and as turns out also here, it is solely sensitive to 
the percolation phenomenon, i.e. to extreme long distance properties. In practice 
the characteristic probability for the occurrence of a network which is nontrivial in 
the three space directions, taking sharp values and 1, allows determination of the 
phases at very low computational cost. 

In addition to showing that the topological properties signal the percolation tran- 
sition it is to be checked wether its transition point coincides with that of the decon- 
fining transition of the fields. This is important because there are examples where 



there is no such coincidence ||18| , |19|, ^]. In the only work |T0[ in which percolation 
is discussed in the context of the deconfinement transition, the coincidence of the 
respective transition points is considered to be an open question. 

It has been observed in Ref. [^ that loops of monopole currents being nontrivial 
in the time direction survive in the deconfinement phase. Such loops have also been 



studied in Ref. [jT4[ by using the concept of a "wrapping number" . This quantity has 



already been introduced in [£T|. It has been shown, because of current conservation, 
to equal the net current flow [|15|. In U(l) gauge theory it is known to fail to 
characterize the transition pT| , |T5|. Thus it appears desirable to investigate in more 
detail what happens in time direction which, as mentioned above, should not be 
involved in the percolation. 

A further question is wether the nontrivial loops in time direction are really simple 
loops and not, in general, networks with a similar number of contacts (crossings of 
monopole currents) as have been observed in U(l) theory. The mean numbers 



per volume of contacts there turn out to be rather insensitive to sizes of lattices 
and of nionopole structures [|15], which indicates a rather uniform density of the 
contacts within the monopole structures. Thus it should be clarified wether the 
same phemomenon occurs also here. 

In the present letter we consider SU(2) gauge theory in the maximally abelian 
projection putting particular effort on careful gauge fixing. We apply the method 
of analyzing networks of monopole currents developed before in U(l) gauge theory 
Ill5| , p!6| and show that the percolation transition of the monopole currents in the 
three space dimensions precisely coincides with the deconfinement transition of the 
fields. In addition we investigate detailed features of the structure of the monopole 
currents in time direction. 

2. Gauge fixing 

The maximally abelian gauge [0 is obtained by performing gauge transformations 
which maximize the quantity 

R = Y.Tria,U,,,asUlJ . (2.1) 

X,fJ, 

The conventional procedure is to perform local gauge transformations iteratively 
throughout the lattice until sufficient accuracy is reached. By introducing an over- 
relaxation parameter u |^ the efficiency of this method can be considerably im- 



proved |2^. For this either the fixed choice a; = 1.7 [^ or choosing stochastically 



^ a; = 1 (no overrelaxation) in 10% of cases and u = 2 otherwise has turned out 
to be advantageous. 

Applying the overrelaxation methods one may still get stuck at some local max- 
imum of (|2.1|) . A method suitable to reach the global maximum is simulated an- 
nealing p^ . We find that for the present purpose to use such a technique is indeed 



crucial. The necessity of using simulated annealing has recently also been pointed 
out in the context of abelian potentials p7[| . 

Our procedure for simulated annealing uses Metropolis sweeps based on a proba- 
bility distribution P{R) ~ exp{aR) in which changes of R by random local gauge 
transformations are proposed. After an appropriate number of sweeps at each a the 
value of a has been increased by a suitable amount and the simulations continued 
at the new value. After the simulated annealing steps in addition an overrelaxation 
algorithm has been applied. Thus essentially first annealing finds the proper maxi- 
mum and then overrelaxation quickly determines the precise numerical value of its 
location. 



To have some control of the quahty of the procedure we have always applied it to 
several gauge copies (created by random gauge transformations of one configuration). 
We then have considered the data separately for the best, the first and the worst 
copy. By varying the quality of the procedure we have checked that it is appropriate 
to associate "best" to the largest value of R and "worst" to the smallest one. 

In standard runs we used 20 a values and 10 sweeps at each of them for annealing. 
We stopped the overelaxation steps after Z = 10~^ has been reached, where Z is 
given by ^ 

z = j-^ui^ir + {x^n , (2.2) 

with X^ = Y.p,{Ui,,xCr:iUl.^ + Ul^^_^,a:iU^,^^-^) , and considered 4 gauge copies. To 
check the reliability of the standard run choices we performed high quality runs with 
600 a-values and 100 sweeps at each of them for annealing. Only minor deviations 
from the best copy results of the standard runs have been observed. These deviations 
have remained within the statistical errors of the measured quantities. 



3. Topological analysis 

In the maximally abelian gauge the coset decomposition of SU(2) with respect 
to the Cartan subgroup U(l) is U^^x = Wfj,,xUfj,^x with tr{Wfi^xO'3) = and where 
Ufj^^x = exp(z^^ ,^.0-3) is diagonal. The abehan physical flux O^i,^^ is given by 9^^^^ = 
0^u,x + '2'^nf,^^x with n^iy^x = 0, ±1, ±2 and 9fj,^^x e [-tt, n) [Q. Defining the monopole 
currents related to the links of the dual lattice by 

'^ ^i,x 7) ^l-iuap\"ap,x+^+u "crp,x+p,) V*^ / 

the conservation law 

Y.iJp,x - Jp,x-p) = (3.2) 

/^ 

has the simple geometrical meaning that incoming and outgoing currents compensate 

at each site of the dual lattice. 

We define current lines in terms of the currents as follows: for J^^^x = there is 
no line on the link, for J^x = =tl there is one line, and for J^x = =t2 there are 
two lines, in positive or negative direction, respectively. Networks of currents are 
connected sets of current lines. For a network N disconnected from the rest the net 
current flow / has the components 

ff^-s = Yl •^M.x for 'Ji^,x e N . (3.3) 



By ( |3.2| ) the net current flows of tlie occurring networks liave to sum up to zero. 

Tlie topological characterization of networks is based on the following observation. 
The elements of the fundamental homotopy group are equivalence classes of paths 
starting and ending at a base point b which can be deformed continuously into 
each other. The generators of the group may be obtained embedding a sufficiently 
dense network N into its space X and performing suitable transformations which 
preserve homotopy. If a given network N does not wrap around in all directions, 
then only the generators of a subgroup are produced. This fact can be utilized for 
an unambiguous characterization of networks. 

In practice we choose one vertex point of N to be the base point b and consider all 
paths which start and end at b. A mapping which shrinks one edge to zero length 
preserves the homotopy of all of these paths. Therefore, by a sequence of such 
mappings we can shift all other vertices to b without changing the group content 
until we finally obtain a bouquet of paths which all start and end at b. In this 
procedure we describe a path by a vector which is the sum of oriented steps along 
the path. The bouquet vectors then form a matrix which is to be analyzed with 
respect to its generator content. This is achieved by a modified Gauss elimination 
procedure which respects current conservation. 



Thus this analysis, on which more details are given in [|T5], |T^, tells in which 
directions the network is nontrivial, which are the directions for which generators 
are found (i.e. for which bouquet loops occur which wrap around the torus). 



4. Results 

Our Monte Carlo simulations have been performed with the Wilson action on 
lattices 8'^ x 4, 12^ x 4 and 18^ x 4. The measurements with gauge fixing and 
analysis of configurations have been separated by 100 simulation sweeps. Typically 
about 100 measurements have been evaluated at each /3-value. 

Figure 1 shows the probability Pnet for the occurrence of a network which is non- 
trivial in the three space directions as function of (3 which we have obtained on 
lattices 18'^ x 4, 12^ x 4 and 8'^ x 4 for best, first and worst gauge copies. The line 
along the best copy results is drawn to guide the eye. It is seen that as soon as the 
lattice is large enough one gets a clear signal, Pnet taking the value 1 in the confining 
phase and in the deconfining phase. In particular, it turns out that for increasing 
lattice size the transition point approaches the one determined from Polyakov loop 
data (where the most recent value is Pc = 2.29895(10) |2^ 



From Figure 1 it can also be seen that going from the data of the best copy via 
the ones of the first copy to those of the worst copy the transition appears to shift 
shghtly to the right. The smallness of this shift shows that the careful gauge fixing 
described in Sect. 2 is successful. Its direction shows that a shift due to bad gauge 
fixing would have the same direction as one due to smaller finite size effects on larger 
lattices. Therefore, reliable gauge fixing is crucial for the present purpose. 

With respect to the number of nontrivial networks we find that in the confining 
phase there is always (in all of our measurements) just one nontrivial network and 
correspondingly its net current flow is zero. In the deconfinement phase there is no 
network which is nontrivial in the three space directions. 

The percolation phenomenon in three dimensions observed here is seen to be com- 
pletely analogous to that found in U(l) gauge theory in four dimensions. Similarly 
as in the U(l) case various alternative possibilities have been checked. It again has 
turned out that precisely the indicated topological properties are the crucial ones 
and give a completely consistent description. 

Next we consider the structure of the monopole currents in time direction. The 
probability for finding a current network which is nontrivial in time direction in the 
P range considered is one on lattices 18^ x 4 and 12^ x 4 , i.e. in these cases we 
find always such a network. On the 8^ x 4 lattice it is one at /9 = 2.2 and decreases 
for larger /3 , at /3 = 2.4 taking a value of about 0.6 for the best gauge copy and of 
about 0.8 for the worst copy. 

Thus it turns out that nontrivial topology in time direction is certain in both 
phases as soon as the lattice is large enough. Again it becomes obvious that the 
8^ X 4 lattice is too small to produce typical results. From the comparison of best 
and worst copies it is again seen that worsening gauge fixing causes a change in 
the same direction as one due to smaller finite size effects on larger lattices. Thus 
careful gauge fixing is important also here. 

A particular feature is that in time direction quite a number of nontrivial networks 
occur in the deconfinement phase. Figure 2 gives the average numbers per lattice 
volume of networks being nontrivial in time direction and of simple nontrivial loops 
in this direction as function of (3 , which we have obtained on lattices 18^ x 4, 12^ x 4 
and 8^ X 4 for best gauge copies. Obviously there is little dependence on the lattice 
size. 

From Figure 2 it is seen that though there is some fraction of simple loops in 
general one has to deal with networks. The observation of a particular fraction at 
a given value of /? is related to the observation already made in U(l) lattice gauge 



theory [|T5| that the average number of contacts of current hues per size has a definite 
value which is rather independent of the sizes of the lattice and of the particular 
current structures. 

The number of contacts at a site is defined by the number of lines arriving at the 
site (or, equivalently, departing from it) minus one. Figure 3a shows the average 
number of contacts of current lines per size ric as function of f3 for all networks and 
Figure 3b Uc for those nontrivial in time direction. It is seen that for the latter ones 
one gets only slightly larger values. 

Thus the situation is analogous as in U(l) gauge theory where the numbers for 



larger networks are only slightly larger than the overall ones |T^. The explanation 
is that the density of contacts gets uniform as soon as the structures are sufficiently 
extended. To demonstrate this effect in the present context we show in Figure 3c 
that a slight increase is also observed if the plaquette-like current loops are omitted. 
It is seen that for the largest lattice then about the same values as in Figure 3b are 
reached. 

The number of contacts turns out to decrease with /3 as one should expect. The 
numerical value at the transition point is about a factor 1.5 smaller as in U(l) theory, 
i.e. the numbers ric are of the same order of magnitude. Thus the observed current 
networks are of the same type. 

If there is more than one nontrivial network, nonzero net current flows become 
possible [0. In the case under consideration they, in fact, occur. Figure 4 shows 
the average numbers per volume of networks with net current flows of modulus 0, 
1 and 2 in time direction as function of P which we have got on lattices 18^ x 4, 
12^ X 4 and 8^ x 4 for the best gauge copy. Apparently |/o| = 1 is most frequent 
while |/o| = 2 is relativly rare. Higher values have not been observed. Again it is 
seen that there is little dependence on the lattice size. 

Possibilities to explain the features occurring in the deconfinement phase by 't 
Hooft-Polyakov monopoles have been discussed in various papers. The situation has 
been reviewed by Smit and van der Sijs [^ who have presented detailed analytical 
investigations. The picture developed by these authors, in fact, predicts nontrivial 
current loops in time direction to persist in the deconfinement phase. Some caution 
with this picture, based on classical solutions and dimensional reduction arguments, 
because of the number of assumptions involved appears, however, appropriate. 

A simpler explanation may focus on the fact that the extension in time direction is 
really very small. Thus, since the current distribution should be sufficiently uniform, 
one can generally expect easy bridging of the torus by monopole currents in that 
direction. 
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Figure captions 

Fig. 1. Probability Pnct as function of P 

on lattices 18^ x 4, 12^ x 4 and 8^ x 4 

for best (circles), first (squares) and worst (diamonds) gauge copies. 

Fig. 2. Numbers per volume Us of networks nontrivial in time direction (circles) 
and of simple loops nontrivial in time direction (squares) as function 
of (3 on lattices 18^ x 4, 12^ x 4 and 8^ x 4 for best gauge copy. 

Fig. 3. Average number per size ric of contacts of current lines, (a) for all networks, 
(b) for networks nontrivial in time direction, (c) for networks except 
plaquet-type ones, on lattices 18^ x 4 (circles), 12^ x 4 (squares) 
and 8^ X 4 (diamonds) for best gauge copy. 

Fig. 4. Numbers per volume Uf of networks with flows |/o| =0 (circles), |/o| = 1 
(squares) and |/o| =2 (diamonds) in time direction as function of (3 
on lattices 18^ x 4, 12^ x 4 and 8^ x 4 for best gauge copy. 
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Figure 2: 
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Figure 3: 
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Figure 4: 
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